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ABSTRACT 

The receptivity to freestream vorticity of the boundary layer 
over a flat plate with an elliptic leading edge is investigated 
numerically. The flow is simulated by solving the incompressible 
Navier-Stokes system in general curvilinear coordinates with the 
vorticity and stream function as dependent variables. A finite- 
difference scheme which is second-order accurate in both space and 
time is used. As a first step, the steady basic-state solution is 
computed. Then a small amplitude vortical disturbance is introduced at 
the upstream boundary and the governing equations are solved time- 
accurately to evaluate the spatial and temporal growth of the 
perturbations leading to instability waves (Tollmien-Schlichting 
waves) inside the boundary layer. Preliminary results for a symmetric, 

2-D disturbance reveal the presence of TS waves aft of the flat- 
plate/ellipse juncture for an aspect ratio 6 leading edge subject to a non- 
dimensional forcing frequency of 230. 


Nomenclature 


U freestream velocity 

L half-thickness of the plate 

v kinematic viscosity 

Re Reynolds number, UL/v 

co frequency of oscillations 

F dimensionless frequency parameter, 

cov/U 2 xl0 6 

AR aspect ratio of the semi-ellipse 
TS Tollmien-Schlichting 


This work is supported by the Air Force Office of Scientific Research and 
NASA Langley Research Center. 
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FS freestream 

X,Y dimensionless Cartesian coordinates parallel and normal to 
the plate, respectively 
V perturbation velocity vector 

u’ streamwise perturbation velocity 

v* normal perturbation velocity 

w’ perturbation vorticity 

A disturbance amplitude coefficient 

T| distance along far-field boundary line 

T| 0 specified mid-point of disturbance 

application line 
o stretching parameter 


1. Introduction 

Boundary-layer stability and transition remains one of the major 
unsolved research areas in fluid mechanics at the time of this 
publication. It is of considerable importance with regard to a myriad of 
technological applications. Linear stability theory has been applied 
extensively in the design of laminar flow control (LFC) systems to 
delay or even eliminate boundary-layer transition. LFC can provide 
major increases in the performance of both aerodynamic and 
hydrodynamic systems (Bushnell and Malik, 1985). Transition 
modeling, which encompasses the accurate prediction of the location of 
transition and the mechanisms which induce it, is crucial to flow 
systems which have both laminar and turbulent regions. In a low- 
disturbance environment, transition in boundary-layer flows is usually 
initiated by the unbounded growth or instability of small disturbances 
in the boundary layer. Boundary-layer stability theory describes the 
development of these self-excited or free oscillations in the laminar 
boundary layer. However, these free or self-excited oscillations are 
usually initiated by some externally forced disturbances such as 
freestream turbulence or sound. 

Receptivity, the process by and through which these external- 
forced disturbances enter the boundary layer and lead to unstable 
waves, thus plays an important role in stability and transition studies. 
Much of what is understood regarding the nature of boundary-layer 
receptivity stems from the high-Reynolds-number asymptotic analysis 
conducted first by Goldstein and colleagues (1983,1985,1989) and 
extended by Kerschen and colleagues(1988,1990,1991). They and their 
colleagues have identified that the conversion of the typically long- 
wavelength freestream disturbances to TS instability waves occurs in 
regions of the boundary layer where there exist rapid local variations 
of geometry and/or local flowfield in the streamwise direction. Much 
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of their work has concentrated on the effect of acoustic disturbances in 
the leading-edge region and regions with short-scale variations in 
surface geometry. Kerschen(1991) has presented results of such an 
analysis for the case of a convected gust. Their analysis has been used 
to aid in the interpretation of experimental results obtained, in general, 
on flat-plate models having elliptic leading edges. Nishioka and 
Morkovin(1986) present a detailed review of some receptivity 
experiments. Recent work by Kendall(1991) in which he studies the 
receptivity problem for the case of freestream turbulence(FST) and 
ongoing work by Parekh et al(1991) also concentrate on similar 
geometries of varying aspect ratio. Kendall has seen evidence of a TS 
structure but experimentalists have, in general, found it difficult to 
isolate a TS component in their measured signal. Much of the difficulty 
may stem from the complex nature of the background flow found in even 
the most sophisticated wind tunnels and therefore in quantifying the 
true disturbance environment (Saric,1990). 

With the advent of powerful supercomputers and numerical 
algorithms, the numerical solution of the full unsteady Navier-Stokes 
equations is fast becoming an integral tool for the study of t e 
receptivity phenomena. One advantage of this technique is that one 
has complete control of the disturbance environment and is therefore 
able to isolate particular kinds of disturbances and study their 
individual effects(Lin,1989). There have been several successful 
numerical works in receptivity to local disturbances, e.g. Fasel et 

In the case of freestream disturbances, Kachanov et al.(19/8) 
solved the incompressible flow over an infinitely thin, semi-infinite 
flat plate, using a finite-difference method to solve the Navier-Stokes 
equations linearized for small disturbances. A von-Karman vortex 
street with its core far from the boundary layer and transverse acoustic 
waves across the leading edge were considered. Murdock (1980) solved 
the full Navier-Stokes equations in parabolic coordinates for flow over 
an infinitely thin flat plate, using a spectral finite-difference method 
and in a 1981 paper presented results for the flow over parabolic bodies 
(with favorable pressure gradient) Gatski and Grosch(1987) solved the 
full incompressible Navier-Stokes equations using finite differences, for 
flow over an infinitely thin, semi-infinite flat plate. Plane 
longitudinal sound waves were considered, however no clear conclusions 
as to the development of the TS wave were presented. Lin (1989) 
numerically simulated the receptivity of the flow over a serm-mfinite 
flat plate with an elliptic leading edge to a plane acoustic disturbance 
for a series of leading edge aspect ratios. His results indicated a 
development of TS waves in the boundary layer, with amplitudes 
which decreased with increasing aspect ratio. Lin et al. (1991) also 
investigated the role of the discontinuity in curvature present at the 
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ellipse flat-plate juncture and found that for a smoothed juncture- 
curvature, the amplitude of the TS waves was roughly halved. 

To our knowledge, unsteady flow over a flat plate including the 
leading-edge curvature has not been solved numerically for the case of 
freestream vorticity. This leading-edge curvature region, with 
inherent pressure gradients and large tangential velocity component, 
has been shown by Lind 989,1 991) to play a vital role in receptivity to 
freestream disturbances. In this work, an attempt is made to 
numerically simulate the receptivity of a two-dimensional boundary 
layer to freestream vorticity and a brief description of the approach 
taken is presented in the next section. 


2. Numerical Formulation 

The present work is an attempt to investigate the initial stage, 
receptivity, of the transition process on a semi-infinite flat plate for an 
incompressible flow, by direct numerical solution of the governing 
Navier-Stokes equations, here cast with vorticity and stream function 
as dependent variables. Following the work of Lin(1989), a C-type 
orthogonal grid is generated around the elliptic leading edge and flat 
plate(Figurel). This mesh is structured in such a way as to allow the 
introduction of a small, localized vortical disturbance at the far-field 
boundary near the basic-state stagnation linefFigures 2a&b), yet 
maintain sufficient resolution of the flat-plate and juncture 
where the development of the instability waves is expected. As the 
surface curvature is changing rapidlyfand is in fact discontinuous at the 
juncture) the mesh is particularly fine in this region. A hyperbolic- 
tangent stretching function is used to insure that no fewer than 12 points 
in the streamwise direction lie within the expected TS wavelength all 
along the flat plate. The equations and boundary conditions are written 
in general curvilinear coordinates and discretized in space and time, 
using second order finite differences. The Modified Strongly Implicit 
Procedure (MSIP) of Schneider and Zedan(1981) is used and is found to 
be robust and stable for reasonably small time steps. 

Initially, the basic-state flow is calculated presuming t e 
following flow conditions: the flow is steady, incompressible with a 
uniform freestream. The basic-state boundary conditions are as follows: 
the no-slip and no-penetration conditions are applied at the wall, a 
symmetric boundary condition along what amounts to the basic-state 
stagnation line (a zero incidence case is considered) is employed, and 
inviscid velocities are specified at the far-field boundaries. At the 
downstream boundary, a homogeneous second-derivative boundary 
condition following Fasel(1976) is used. 
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The basic-state, upon application of the afore-mentioned 
boundary conditions, is solved using a transient approach and spatially 
varying time steps. Once a basic-state solution is obtained, a small- 
amplitude (u' ma x/U=4xlO- 3 ), periodic, localized disturbance is 
applied as an unsteady boundary condition at the far-field boundary. 
The resulting unsteady flow and the temporal and spatial development 
of the perturbations are determined by solving the governing equations 
subject to this modified boundary condition time accurately. Current 
results are presented in the following section. 


3. Results and Discussion 

The steady-state solution is obtained for an aspect-ratio(AR)=6 
semi-ellipse. The Reynolds number(Re), based upon the half-thickness 
of the leading edge and the freestream velocity U is 2400. The 
downstream boundary of the computational domain is chosen such that 
the branch I of the neutral stability curve (according to linear parallel 
theory) lies well within the computational domain. The farfield 
computational boundary is placed such that it is 20 boundary-layer 
thicknesses from the plate in the downstream region. The mesh size 
used in the present analysis is 260 points in the tangential direction and 
100 points in the normal direction, although a calculation for a mesh of 
360x100 yielded an identical result for the basic-state calculation. 
Wall vorticity is presented in Figure 3. Nearly Blasius flow is 
recovered in the flat-plate region as evidenced by the streamwise 
velocity profiles depicted in Figure 4. 

The unsteady calculation was carried out for a dimensionless 
frequency parameter F=230. The forcing amplitude applied has a 
maximum streamwise velocity perturbation of (u’max/U)=4xlO • 
disturbance was introduced over a small region just above the stagnation 
line along the farfield boundary in the form of a periodic boundary 
condition cast in the following form : 



4(7?- Tip) -2 
< 7 


~(v~Vo y t** 


sin(£or) 


( 1 ) 


Note that this represents the vorticity associated with the 
introduction of a perturbation of velocity normal to the farfield 
boundary. The velocity perturbation is specified such that it is 
symmetrical about a specified location along the farfield boundary, 
rio , so that no net mass is introduced into the domain. The disturbance 
was found to propagate downstream, as expected, at the convective 
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speed and as such the wavefront of the disturbance reached the 
downstream boundary during the 6th cycle of forcing. Figure 5 is a 
presentation of instantaneous disturbance streamlines observed at the 
end of a series of cycles of forcing which reveal the quasi-steady 
structure of the solution over much of the plate following 6 cycles. A 
temporal Fast Fourier Transform (FFT) analysis conducted over several 
cycles also confirmed the quasi-steady nature of the flow after 6 cycles. 
In figure 6 typical u'-velocity profiles are displayed and in figure 7 the 
change of u' with X is presented. From these figures it is evident that 
the signal is a complex combination of more than one wave. In an 
attempt to extract some information regarding the specific components 
of this signal, some educated assumptions were applied. Kerschen 
points out that the freestream vorticity will decay exponentially as 
one approaches the wall. This is evident in the streamwise u 
plotsffigure 7). Note the change in signal wavelength with distance 
from the wall. Near the wall, the signal wavelength is less than half 
of that found in the freestream and, qualitatively, the wavelength of 
the signal appears to approach that of the TS wave. Thus, to better 
assess the nature of the disturbance, spatial correlations with signals of 
freestream and predicted TS wavelength were conducted. In figure 8, a 
series of correlations between the total signal and a signal with a 
wavelength based on forcing frequency and freestream velocity at 
different Y-locations is presented. Clearly, this correlation supports 
our supposition that the direct contribution of the freestream 
perturbation to the total signal near the wall is minimal. If one then, 
as a first guess, presumes the signal wavelength one point from the wall 
to be that of the TS wavefa length as it turns out very close to that 
which is predicted by linear theory), then one can use the correlation 
with this wavelength to determine the phase of this portion of the 
signal. Phase information for this portion, hereafter referred to as the 
TS portion of the total signal, is presented in figure 9 for a 
representative X-location aft of the juncture Note the 180° phase shift 
which is associated with the TS wave. A comparison with the total 
signal phasefobtained from a temporal FFT over one cycle of forcing) 
with that obtained for the freestream and TS components is displayed 
in figure 10. Note that the total signal and the TS’ are in-phase near 
the wall, as expected, but as one moves away from the wall, the 
freestream portion of the signal begins to dominate. The correlation can 
also give a rough estimate of the signal amplitude, which again 
displays the characteristic TS mode shape (Figure 11). If one uses the 
phase information obtained from the spatial correlations in concert 
with the temporal FFT of the total signal, an estimate of the 
respective component amplitudes which accounts for signal modulation 
with X may be obtained. Amplitude plots are presented in Figure 12 for 
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a series of X-stations using this technique. As a measure of the 
receptivity of the disturbance, one may compare the ratio of the 
maximum TS amplitude to that of the disturbance in the freestream. 
From this analysis we find that for an AR=6 leading edge at a forcing 
frequency F=230, the receptivity is of order 10 


4 Conclusions 

Preliminary results from a numerical simulation of the 
receptivity of a two-dimensional laminar boundary layer to a 
symmetric vorticity wave in the freestream suggest the presence of TS 
waves aft of the ellipse-flat-plate juncture. While the signal 
decomposition technique is imprecise, the receptivity, based on 
maximum TS amplitude to freestream amplitude, is estimated to be of 
order 10' 1 . In addition to a more detailed analysis of the current result 
to include a comprehensive grid study and an analysis of the leading 
edge flowfield, studies of the effect of disturbance position, amplitude, 
frequency, leading-edge aspect ratio, and the effect of juncture 
curvature are planned. Similar cases will be investigated for the 
particularly interesting case of an asymmetric disturbance. Also in 
progress is a repeat of this case on a computational domain which 
extends past branch II. 
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Figure 1. C-grid used in conputations. 
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Figure 2a. Blow-up of disturbance input region. 
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Figure 2b. Schematic of disturbance velocity profile at input location. 
(V has both a u’ and v' component). 



Figure 3. Basic state wall vorticity. 
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Figure 5. Disturbance 
after 5th, 6th anc 7th 


Streanrfunction contours 
cycles of forcing. 
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Figure 7. Streamwise-perturbation velocity profiles 
after 9 cycles of forcing. 
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Figure 8. Spatial correlation of the total signal with a signal of 
freestream wavelength. 
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5. RESOURCES AND PERSONNEL 

One of the principal strengths of our team at Arizona State University is its broad skills in 
analysis, computations, and experiments. We facilitate day-to-day communication between the 
computational work and the experimental work through two IRIS Graphics Workstations (3030 
and 3130) and two DEC 5000 Workstations. The system, with state-of-the-art, real-time, three- 
dimensional, color-graphics software (PLOT3D), is equipped with an extensive multi-user and 
multi-task environment with twelve serial lines. Users are able to share the same data base or 
experimental information. This provides the heart of the interaction of the analytical, 
computational, and experimental research. 

In addition to the super computers at NASA facilities and Princeton/NSF Consortium, the 
network includes access to the IBM 434 1/VM and Harris/VS computers, the IBM 3090 Class VI 
machine, and the Cray XMP on campus as well as the MASSCOMP. The College of 
Engineering at ASU is currently also equipped with several VAX/780 and VAX/785 
minicomputers exclusively for research purposes (each office and laboratory has a hard-wired 
RS232 interface). These minicomputers are excellent systems for program development. The 
IRIS and DEC machines can access all the features available in those minicomputers through the 
existing local area networking (Ethernet) on the campus. Furthermore, the system can 
communicate directly with NASA research facilities to share information through telephone 
couplings. The full array of computer capabilities from super-mini to super-super is in place for 

the research. 

Ray-Sing Lin, a PhD student, and Nay Lin, a PhD student, are the two candidates for the 
Fellowship. Ray, Nay, and Tom Buter (also a PhD student) each spent one month of the summer 
of 1991 at Langley Research Center working on these and similar problems, interacting with 
NASA and ICASE personnel. 

The principal investigator for this work is Helen L. Reed, Associate Professor of 
Mechanical and Aerospace Engineering. Professor Reed has spent the last ten years conducting 
theoretical and computational research on problems of boundary-layer stability specifically 
applied to the ACEE/LFC programs. She also spent one month of the summer of 1991 at 
NASA/Langley Research Center. Her resume is attached as Appendix I. 
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must take in the preparation and training of the U.S.'s future scientists and lay people. 

Member of National Academy of SciencesiN ational Research Council Aerodynamics Panel which is a pari of 
the Committee on Aeronautical Technologies of the Aeronautics and Space Engineering Board , 
Commission on Engineering and Technical Systems , November 1990-March 1992. This is the advisory 
group to NASA and the U.S. Congress concerning the directions NASA must take in order to enable the 
U.S. to remain competitive in the world arena. 

Member of U.S. National Transition Study Group under the direction of Eli Reshotko, 1984-Present. 

Associate Editor, Annual Review of Fluid Mechanics , 1986-Present. With the work divided equally among 
John Lumley, Milton Van Dyke, and myself, we are responsible for complete selection, refereeing, and 
editorial correction of all articles in each yearly volume. 

Originator of Annual Picture Gallery of Fluid Motions at annual meetings of the American Physical Society , 
November 1983, Houston. Responsibility for the gallery at American Physical Society meeting at Brown 
University, November 1984: Eugene, November 1987: Buffalo, November 1988, Palo Alto, November 
1989: Cornell, November 1990: Scottsdale 1991. 

Member of AlAA Technical Committee on Fluid Dynamics , 1984-1989. APS/AIAA and ASME-AMD/AIAA 
Liaison. Prepared Aerospace Audit article for Fluid Dynamics, Aerospace America , December 1988 (did 
not appear because of difficulties between Technical Committee and editors). 

Member of Fluid Mechanics Technical Committee of the Applied Mechanics Division of the ASME , 1984 - 
Present. 

Member of Steering Committee for National Fluid Dynamics Congress , June 1988-Present. 

University Representative of National Science Foundation Fluids Engineering Workshop. Group Leader and 
Coordinator of final document for Unsteady Flow Panel, Savannah, Georgia, September 17-20, 1986. 


H.L. Reed: Swept-Wing Flows (NASA-Langley, Graduate Program in Aeronautics) 
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Chairperson of ASME Junior A war ds Committee 1989-Prescnt. Vice-Chairperson 1987-1988. 

Member of Scientific Committee of 1992-3 IUTAM Symposium on Nonlinear Stability of Nonparallel Flows, 


January 1990-Present. c ,qo/ 

Chairperson, 2nd Annual Arizona Fluid Mechanics Conference, Arizona State University, April 4-5 1986. 
Technical Chairperson, AJAA 19th Fluid, Plasma Dynamics, and Lasers Conference , ; Honolulu June 9 > • 

Co-Chairperson, with Dan Jankowski, 44th Annual American Physical Soaety/Dmsion of Fluid Dynamics 


Meeting, Scottsdale, November 1991. n r _ 

Co-Organizer with Dr. Daniel Reda (Sandia National Laboratones) of Symposium 
Theoretical/Numerical Studies of Boundary-Layer Stability and Transition 
ASME/JSME Fluids Engineering Conference, Portland, June 23-26, 1991. 

Chair of 21 various symposia and sessions at international conferences. 


on Experimental and 
at the First Joint 


CONSULTING 

Westinghouse Electric and Naval Underwater Sea Center (1981) 
Pratt and Whitney (1986-1991) 

ICASE Consultant, NASA/Langley Research Center (Current) 
Colorado Research Development Corporation (Current) 



